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Abstract 

We derive exact renormalization-group recursion relations for an Ising 
model, in the presence of external fields, with ferromagnetic nearest- 
neighbor interactions on Migdal-Kadanoff hierarchical lattices. We con- 
sider layered distributions of aperiodic exchange interactions, according 
to a class of two-letter substitutional sequences. For irrelevant geometric 
fluctuations, the recursion relations in parameter space display a nontriv- 
ial uniform fixed point of hyperbolic character that governs the universal 
critical behavior. For relevant fluctuations, in agreement with previous 
work, this fixed point becomes fully unstable, and there appears a two- 
cycle attractor associated with a new critical universality class. 

PACS number(s): 05.50; 64.60; 75.10.H 

The discovery of quasi- crystals has motivated the investigation of the effects 
of geometric fluctuations, associated with aperiodic interactions, on the critical 
behavior of spin systems JD. According to the heuristic Harris criterion it is 
known that fluctuations produced by disordered interactions change the critical 
behavior of ferromagnetic systems whose uniform counterparts are characterized 
by a diverging specific heat at the transition. Although aperiodic interactions 
come from deterministic rules, and are thus highly correlated, they also intro- 
duce extra fluctuations that may lead to drastic changes in the critical behavior. 
As proposed by Luck[||, there is indeed a criterion of relevance of the geometric 
fluctuations, which has been confirmed by calculations for quantum spin chains 
with aperiodic exchange interactions Q. 

Taking advantage of the structure of the Migdal-Kadanoff hierarchical lat- 
tices, it is relatively easy to write exact recursion relations to analyze the critical 
behavior of ferromagnetic Ising and Potts models, in zero field, with layered 
aperiodic exchange interactions that are chosen according to a large class of 
substitutional sequences (|, @, ||. There is an exact analog of the Harris-Luck 
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criterion to gauge the relevance of the geometric fluctuations in terms of the 
wandering exponent of the substitutional sequence, the number of states (of 
the Potts model), and some structural features of the particular hierarchical 
lattice|| For irrelevant fluctuations, the (universal) critical behavior is gov- 
erned by a nontrivial uniform fixed point in parameter space of the recursion 
relations. For relevant fluctuations, at least in the case of two-letter substitu- 
tional sequences, this nontrivial uniform fixed point becomes fully unstable, and 
there appears a two-cycle hyperbolic attractor that gives rise to a new class of 
universal critical behavior |||. 

In the present publication, we report an extension of these calculations to 
analyze the field behavior of an Ising model with layered aperiodic exchange 
interactions on Migdal-Kadanoff hierarchical lattices. In one dimension, we 
regain results of previous work by Achiam, Lubensky, and Marshall || for the 
ferromagnetic Ising chain with aperiodic exchange interactions (according to 
the standard two- letter Fibonacci sequence). In the uniform limit, we recover 
the results of Bleher and Zalys(lCj for the field behavior of the Ising model 
on the diamond hierarchical lattice. It should be mentioned that spin models 
with uniform interactions on general hierarchical lattices have been thoroughly 
investigated in a series of papers by Kaufmann and Griffiths Jli[. 

As in the work of Pinho, Haddad, and Salinas|| Ej, we initially consider a 
period-doubling two-letter sequence, 

A AB, (1) 

B -> AA, (2) 

so that the successive application of this inflation rule on letter A produces the 
sequence 

A — > AB — » ABAA -> ABAAABABA -+ ... (3) 
The Hamiltonian of the nearest-neighbor Ising model is written in the form 

n = - ^2 [JijWj + \h%A a i + °j) + \ h h( a i - a o)) ' ( 4 ) 

where a. L = ±1 is a spin variable on site i of a hierarchical diamond structure, 
{ Jij} is a set of exchange interactions, hij and hh are external fields associated 
with the bonds of the structure, and the sum is over all pairs of nearest- 
neighbor sites. As in the one-dimensional problem ||, it is essential to take into 
account the staggered external field h' in order to avoid any inconsistency in 
the final recursion relations. In Fig. 1, we show three successive stages of the 
construction of a diamond lattice (bond length 6 = 2, and number of branches 
m = 2) for the period-doubling sequence given by Eq. (0). As indicated in the 
figure, we construct a layered distribution of both the exchange parameters, J a, 
and Jb, and the external fields, Ha, h/ A , hs, and h B , along the branches of the 
simple diamond lattice. By decimating the intermediate spin variables (see Fig. 



2 



1), it is straightforward to write the recursion relations 

/ _ XAXB{l + y 2 AyB z A z B) + VAVBZAZb(x 2 b + x\) 

A x A x B {l + y 2 A y 2 B z 2 A z 2 B ) + vaVbzazb^ + x 2 B x 2 A ) ' 



(5) 



\{l + y A z\) + y\z A {l + x A y 



(6) 



, _ VAyB{x A x B + vaVbzazb) /-s 

y A 77~, \) {') 

z A z B (l + xaxb yA yB z a z b ) 



I _ VA^A + V\ Z A) (a\ 
UB „2 H I T 2 ,,2 7 2 \ ' \°) 

iU + x aVa / 



VaZb{x b + XAyAVBZAZB) 
VBZa(xa + X B yAyBZAZB) ' 



(9) 



and 



4 = 1, (10) 



where ia,b = exp (-2(3Ja,b), Va,b = exp (-/3/m,b), = exp 
and /3 is the inverse of the temperature. 



A,B 



In the uniform case, Ja = Jb = J, h-A = h,B = h, and /i^ = /i^ = /i T , Eqs. 



(§)-© reduce to the recursion relations 

, a; 2 (l + y 4 z 4 ) + 2a; 2 y 2 z 2 



x 2 (l + y 4 z 4 ) +y 2 z 2 (l +x 4 )' 



(11) 



and 



, _ y 2 x 2 + y 2 z 2 , . 

^ z 2 l + x 2 y¥' 1 j 



(13) 



which can be compared with the expressions obtained by Bleher and ZalyspOf 
for the simple Ising model in a (uniform) field. For the aperiodic model, but in 
the absence of the external fields, we recover the equations obtained by Pinho et 
al.|| ^j. Due to the presence of a staggered field, it is important to see that the 
form of the recursion relations depends on the way we label the bonds on the 
diamond lattice. To be consistent with the calculations of Bleher and Zalys, we 
assume a positive staggered field at the intermediate spins (<Xj) of the polygon, 
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and an opposite field at the nodes (<7j). With a change of order, there is no 
automatic cancellation of the staggered field associated with the B bonds. 

The nontrivial uniform fixed point of the recursion relations (0)- ( |To| ) is 
given by yA = l/B = za = zb = corresponding to an Ising model in zero field 
(fiA = hi} = h A = hg = 0), and 

x* A = x* B = x = 0.295598..., (14) 

which comes from the cubic equation x 3 + x 2 + 3x — 1 = 0, and is the usual result 
for the uniform model. The trivial fixed points are located at yA = Vb = za — 
Zb = 1, with x* A — x* B = or 1. The linearization of the recursion relations 
about the nontrivial uniform fixed point, yields the block-diagonal form 
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(15) 



where C = (1 — x 2 )/(l + x 2 ) = 0.839286.... The thermal eigenvalues, given by 

Ai = 2C = 1.67857..., and A 2 = -C = -0.839286..., (16) 

are identical to the findings of Pinho, Haddad, and Salinas||, in zero external 
field. The magnetic eigenvalues are given by 

A 3 = 2 (1 + C) = 3.678572..., A 4 = -2, (17) 



A 5 = -C = -0.839286..., and A 6 = 0. (18) 

For the Ising model on the diamond lattice with a uniform distribution of ex- 
change interactions (as in the work of Blchcr and Zalys), the linearization about 
the nontrivial fixed point of Eqs. (|ll])-([l^) leads to the set of three uniform 
eigenvalues, Ai = 2C = 1.67857..., A 3 = 2 (1 + C) = 3.678572..., and A 6 = 0. 
Therefore, as in the work of Pinho, Haddad, and Salinas, and in agreement 
with the analog of the Harris-Luck criterion, we sec that the geometric fluctu- 
ations associated with the period-doubling sequence are unable to change the 
ferromagnetic critical behavior of the Ising model on the diamond lattice. The 
absolute value of the smaller thermal eigenvalue is less than 1. The absolute 
values of the eigenvalues associated with the uniform field are both larger than 
one. The staggered field is irrelevant as in the uniform case. Also, it should be 
noted that the same kind of behavior still holds for a more general diamond-type 
lattice, with 6=2 bonds and m > 2 branches. 

Due to the peculiar structure of the diamond hierarchical lattice, there is 
a change in the magnetic eigenvalues if we make another choice for the sign of 
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the staggered field. This effect is already present in the model with uniformly 
distributed exchange interactions. In fact, if we linearize the new recursion 
relations of the uniform model, in direct and staggered applied fields, we obtain 
the same values for Ai and A3, but Ag = —2, instead of A$ = 0. For the 
aperiodic model, we obtain the new values A4 = —1.78151..., A 5 = —0.942221..., 
and A 6 = —2. Hence, we see that we cannot attribute this change of magnetic 
eigenvalues to the introduction of aperiodic exchange interactions. It is evident 
that aperiodicity is not the cause for these changes. 

To give an example of relevant geometric fluctuations, we now consider the 
period-three two-letter sequence, 

A ABB, (19) 

B AAA. (20) 

The nearest-neighbor Ising model is defined on a Migdal-Kadanoff lattice with 
6 = 3 bonds and m = 2 branches. Although somewhat tedious, it is straightfor- 
ward to write the recursion relations 



/ _ Fi (xA,XB,yA,UB, ZAZ B ) F 2 (xA,XB,yA,yB,ZA,ZB) 
JA F 3 {XA,XB,yA,yB,ZA,ZB)F 4 (xA,XB,yA,yB,ZA,ZB) 



(21) 



/ _ Fi (x A , xa, VA, Va, z a za) F 2 (x a , x a , Va, Va, zaza) < 22 s 

B Fs(XA,XA,yA,yA,ZA,ZA)F 4 (XA,XA,yA,yA,ZA,ZA)' 



1 _ yAVBZAFj [x A , Xb 1 yA ; Vb ; Z A . Zb) , s 

Z B F 3 {XA,XB,yA,yB,ZA,Z B ) 

, _ y A F 4(xA,XA,yA,yA,ZA,ZA) , OA , 

F 3 (x A , x A , va, va, zaza) 



, yAZAZBFi(x A ,x B ,yA,yB, za.zb) 

Z a — T7~i \ , \^) 

yBF 2 (x A ,XB,yA,yB,ZAZB) 



, z\Fi{x A ,X A ,yA,yA,ZAZA) 

Z B — '1 \ — ) \^°) 

F 2 {XA,X A ,yA,yA,ZA,ZA) 



where 



Fi (x A ,XB,VA,yB,ZAZB) = X A ZAZB + XaXbUbZAZb +X B yAVB + X B yAy B z B, 

(27) 

F 2 (x A ,XB,VA,yB,ZAZB) = XbZaZb + x A x B y A y B + x B y B z A z B +x A y A y%z B , 

(28) 
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F 3 (XA, XBiUA, VB,ZA,Zb) = Z A Z B + XAXBVAVB + X A X B yAy\z 2 B + X 2 B y 2 B ZAZ%, 

(29) 

and 

F 4 (xa ^Bj Va ) Vb j z a .z b ) = x a x b zazb + x A x B y 2 B z A z B + x B y A yB + VaV%z%. 

(30) 

Besides the trivial fixed points, there is a nontrivial (uniform) fixed point, 
given by y a — Vb — za = z B — 1 (zero direct and staggered fields), and 



x\ = X * B = x = 0.119726..., 



which comes from the equation 



l + 3x 2 



(31) 



(32) 



and is identical to the fixed point associated with the uniform model. We now 
linearize the recursion relations about this nontrivial fixed point. Again, we have 
a block-diagonal form involving separate thermal and magnetic contributions, 
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F 



3 + x 2 



l + 3x 2 

We then have the thermal eigenvalues, 

Ai = 3D = 1.854101..., and A 2 = -2D = -1.236067. 



(35) 
(36) 
(37) 



which are identical to the findings of Pinho, Haddad, and SalinasQ, in zero 
external field, and the magnetic eigenvalues, 



A, 



2E 
~F 



5.780107... 



A 4 



-A(E + x 2 F) = -3.854101..., 



(38) 
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A 5 = 0, and A 6 = — = 0.692028.... (39) 
E 

For the Ising model on the diamond lattice with a uniform distribution of ex- 
change interactions (as in the work of Bleher and Zalys), the linearization about 
the nontrivial fixed point of Eqs. (^l|)-(|2^) leads to the set of three uniform 
eigenvalues, Ai = 3D > 1 , A 3 = 2E/F > 1, and < A 6 = 2F/E < 1 (the 
external staggered field is irrelevant). As in the work of Pinho, Haddad, and 
Salinas Q, since Ai > 1 and |A 2 | > 1, the uniform fixed point is unstable along 
both thermal directions, and cannot be reached from physically acceptable ini- 
tial conditions. This type of behavior, in full agreement with the exact analog 
of the Harris-Luck criterion, still holds for more general lattices with b — 3 but 
m > 2. 

In the more recent analyses of Haddad, Pinho, and Salinas^, in addition to 
the (unstable) uniform fixed point, there is a two-cycle attractor that governs 
the critical behavior. Recursion relations (|2l|)-(|2^) are also associated with a 
two-cycle attractor, located at tja — ys = za = zb = 1, and 

x\ = 0.155117..., x* B = 0.00740154..., (40) 

and 

x* A = 0.0287405..., x* B = 0.191409.... (41) 

The linearization of the second iterates of the recursion relations about this 
two-cycle attractor leads to the thermal eigenvalues, 

Ai = 3.255710..., and A 2 = 0.311746..., (42) 

where | A.2 1 < 1 indicates that the attractor can be reached from physically 
acceptable initial conditions, and the magnetic eigenvalues 

A 3 = 34.384598..., A 4 = 15.008653..., (43) 

and 

A 5 = 0.267660..., and A 6 = 0, (44) 

which indicates the irrelevance of the applied staggered field. These results are 
consistent with the recent work of Haddad, Pinho, and Salinas ||. The nontriv- 
ial diagonal fixed point cannot be reached. In zero field, a new universal critical 
behavior is characterized by the eigenvalues about the cycle-two attractor. 

In conclusion, we have investigated a ferromagnetic Ising model with ape- 
riodic exchange interactions, in the presence of direct and staggered external 
fields, on a class of Migdal-Kadanoff hierarchical structures. The geometric 
fluctuations introduced by the aperiodic interactions may change the character 
of the fixed point associated with the critical behavior. For relevant geometric 
fluctuations, the critical behavior belongs to a new universality class associated 
with a two-cycle attractor. The inclusion of a magnetic field does confirm these 
findings. 

The authors thank Fundacao de Amparo a Pesquisa do Estado de Sao Paulo 
for financial support. 
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Figure 1: Some generations of the simple diamond hierarchical lattice (6 = 2 
bonds and m = 2 branches). The layered exchange interactions are chosen 
according to the period-doubling two-letter substitution sequence, A — > AB, 
B AA. 
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